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ABSTRACT 

We report molecular line and dust continuum observations, made with the 
SEST telescope, towards four young high-mass star forming regions associated 
with highly luminous (£ > 6 x IO^Cq) IRAS sources (15290-5546, 15502-5302, 
15567-5236 and 16060-5146). Molecular emission was mapped in three lines of 
CS (J=2^1, 3^2 and 5^4), two lines of SiO (J=2^1 and 3^2), two rotational 
transitions of CH3OH (Jfc=3fc 2k and 2^ —>■ 1^), and in the C^''S(J=3— >2) line. 
In addition, single spectra at the peak position were taken in the C0(J=1— >0), 
i^CO(J=1^0) and C^«O(J=1^0) hues. We find that the luminous star form- 
ing regions are associated with molecular gas and dust structures with radii of 
typically 0.5 pc, masses of ~ 5 x 10^ M©, column densities of ~ 5 x 10^^ cm~^, 
molecular hydrogen densities of typically ~ 2 x 10^ cm~^ and dust temperatures 
of ~ 40 K. The 1.2 mm dust continuum observations further indicate that the 
cores are centrally condensed, having radial density profiles with power-law in- 
dices in the range 1.6 — 1.9. We find that under these conditions dynamical 
friction by the gas plays an important role in the migration of high-mass stars 
towards the central core region, providing an explanation for the observed stellar 
mass segregation within the cores. 

The CS profiles show two distinct emission components: a bright component, 
with line widths of typically 5 km s~^ (FWHM), and a weaker and wider ve- 
locity component, which extends up to typically ±13 km s^^ from the ambient 
cloud velocity. The SiO profiles also show emission from both components, but 
the intensity of the pedestal feature relative to that of the bright component is 
stronger than for CS. The narrow SiO component is likely to trace warm ambient 
gas close to the recently formed massive stars, whereas the high velocity emis- 
sion indicates mass outflows produced by either the expansion of the HII regions. 
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stellar winds, and/or coUimated outflows. We find that the abundances of CS, 
CH3OH and SiO, relative to H2, in the warm ambient gas of the massive cores 
are typically 4 x 10~^, 6 x 10~^, and 5 x 10~^^, respectively. 

Subject headings: ISM: clouds — ISM: dust — stars: formation — stars: massive 
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1. INTRODUCTION 

During the last decades several observational efforts have been carried out to determine 
the characteristics of the molecular gas and dust associated with massive star forming 
regions. Surveys of molecular emission in high density tracers (Plume et al. 1992; Juvela 
1996; Plume et al. 1997; Shirley et al. 2003) and of dust continuum emission (Beuther et 
al. 2002; Mueller et al. 2002; Faiindez et al. 2004; WiUiams et al. 2004; Garay et al. 2007), 
made with single dish telescopes, show that the maternities of high-mass stars are regions 
of molecular gas with distinctive physical parameters, with typical radii of 0.4 pc, densities 
of 5 X 10^ cm~'^ and masses of 4 x 10^ Mq. The dust observations also show that these 
regions have dust temperatures of typically 32 K, indicating that a luminous energy source 
has been already formed inside them. These structures, also refereed as massive and dense 
cores, have typically column densities of 3 x 10^'^ cm~^ (Garay et al. 2007) which make 
them very dark at optical wavelengths and even dark at infrared wavelengths (Menten et 
al. 2005, Rathborne et al. 2006). The IRAS luminosity of most regions in these studies are 
in the range 1 x 10"^ — 4 x 10^ Lq. 

In order to investigate possible differences, in the physical and chemical characteristics 
as well as in the stellar population, between massive and dense cores of different luminosities, 
we have carried out a multi-wavelength study of four highly luminous IRAS sources in 
the southern hemisphere selected from the list of IRAS point sources with colors of UC 
HII regions and detected in CS(2— i>l) emission compiled by Bronfman et al. (1996). The 
observed objects are listed in Tabled! Cols. 2, 3, and 4 give, respectively, the name, right 
ascension, and declination of the associated IRAS source. The kinematic distance to the 
objects is given in col. 5. The four objects are located in the tangent region of the Norma 
spiral arm and therefore there is little ambiguity in the kinematic determination of their 
distances; their vlsr^ are nearly the maxima allowed by pure circular motion in their lines 
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of sight. The total far infrared luminosities, Lpm, are given in col.6. The objects have Lpm 
greater than 6 x 10^ Lq and their CS(2^1) line profiles show broad molecular line wing 
emission suggesting the presence of powerful molecular outflows. 

Here we report molecular line and dust continuum observations towards these four 
regions. The aim of these observations is to undertake a detailed study of the molecular 
gas and dust content of luminous massive and dense cores in order to determine their 
physical conditions, kinematics, and chemistry. Molecular emission was mapped in a 
number of species that are commonly used as diagnostics of the conditions of massive and 
dense cores: carbon monosulfide, silicon monoxide, and methanol (e.g., Zinchenko et al. 
1995; Juvela 1996). Thermal emission of SiO is both a good probe of high velocity gas in 
regions of massive star formation (Downes et al. 1982) and a tracer of high-temperature 
chemistry (Ziurys et al. 1989). The methanol abundance is highly enhanced in the hot 
gas surrounding infrared sources (Menten et al. 1988) as well as in shocked gas (Liechti & 
Walmsley 1997), thus thermal methanol emission is an excellent probe of the environments 
of high-mass YSO's. The CS and dust continuum observations are intended to determine 
the density structure of the cores. In a subsequent paper (Chavarria et al. 2009; Paper II) 
we will present the characteristics of the stellar content within these regions. 



2. OBSERVATIONS 



The molecular line and dust continuum observations were carried out using the 15-m 
Swedish-ESO Submillimetre Telescope (SEST) located on La Silla, Chile. 
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2.1. Molecular lines 

The molecules, transitions, and frequencies observed and the instrumental parameters 
are summarized in Table [2l 

The CS(2 1) observations were made between December 1990 and January 1992. 
During this period the telescope was equipped with a cooled Schottky mixer and an 
acousto-optical spectrometer with 2000x43 KHz channels. The typical system temperature 
above the atmosphere was 500 K. The CS(2 1) line emission was mapped with spacings 
of 45" within regions of at least 135" x 135". The observations were done in the frequency 
switch mode with a frequency throw of 15 MHz. 

The SiO(2 ^ 1), SiO(3 ^ 2), CH30H(2fe ^ U), CH30H(3fc 2^), CS(3 ^ 2) and 
C^^S(3 — s> 2) observations were made during March of 1997. In this epoch the telescope was 
equipped with SIS receivers, and it was possible to simultaneously observe lines at 2 and 
3 mm wavelengths. Single-sideband receiver temperatures were typically 120 K for both 
receivers. As backend we used high resolution acousto-optical spectrometers providing a 
channel separation of 43 kHz and a total bandwidth of 43 MHz. The emission in the above 
lines was mapped within regions of at least 90" x 90", with 30" spacings. The observations 
were performed in dual beam- switching mode, with a beam separation of 11'47" in azimuth. 
The integration times on source per position and the resulting rms noise in antenna 
temperature are given in Table [3l Within the available bandwidths, three rotational lines 
of CH3OH could be observed at 2 mm (3o 2o A+, 3_i 2_i E, and 3o ^ 2o E lines) 
and four rotational lines at 3 mm (2i li E, 2q —>■ 1q E, 2q —>■ 1q A'^, and 2_i —>■ l_i E 
lines). During March 2001 we further observed the SiO(5 4) (217104.935 GHz) and 
CH30H(5_i 4o) (84521.21 GHz) hues toward the center of the G330.949-0.174 core. 
The receivers and backends were the same as described above. Single-sideband receiver 
temperatures were 1200 K at 217 GHz and 190 K at 85 GHz. 
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The CS(5 4) observations of G324.201, G329.337, and G330.949 were made between 
September 1991 and January 1992. In this epoch the telescope was equipped with a cooled 
Schottky mixer, and an acousto-optical spectrometer with 2000x43 KHz channels; however, 
a setup of 1000x86 KHz resolution elements was used. The typical system temperature 
above the atmosphere was 1850 K. The CS(5 4) line emission was mapped with spacings 
of 15" within regions of at least 60"x45". The observations were done in the frequency 
switch mode with a frequency throw of 15 MHz. The CS(5 4) observations of G328.307 
were made in May 2001, with an SIS receiver and the acousto-optical spectrometer with 
2000x43 KHz channels. Typical system temperature above the atmosphere was of 450 K. 
Observations were done in dual beam switching mode, with a beam separation of 11' 47" in 
azimuth, in a grid of 9 positions with spacings of 45". 

Observations of the CO, ^^CO, and C^^O emission in the J=l— >0 lines were made 
during May, 2000 and June, 2001. Spectra were obtained only toward the peak position 
of the CS emission in each region. These observations were performed in the position 
switched mode, with OFF positions taken from Bronfman et al. (1989). The goal of these 
observations was to obtain the data needed to determine CO column densities, which in 
turn will permit to derive abundance enhancements relative to CO. 



2.2. Dust continuum 

The 1.2-mm continuum observations were made using the 37-channel SEST Imaging 
Bolometer Array (SIMBA) during October, 2001 and August, 2003. The passband of the 
bolometers has an equivalent width of 90 GHz and is centered at 250 GHz. The HPBW of 
a single element is 24" and the separation between elements on the sky is 44". We observed 
in the fast mapping mode, using a scan speed of 80" s"-*^. Each observing block consisted 
of 50 scan lines in azimuth of length 800" and separated in elevation by 8", giving a map 
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size of 400" in elevation. This block required ~ 15 minutes of observing time. We observed 
2 — 4 blocks per source, achieving rms noise levels of typically 50 mJy beam~^. The data 
were reduced according to a standard procedure using the software package MOPSl, which 
included baseline subtraction and rejection of correlated sky-noise. Flux calibration was 
performed using a sky-opacity correction and a counts-to-fiux conversion factor derived 
from maps of Uranus. Uncertainties in the absolute calibration and pointing accuracy are 
estimated at 20% and 3", respectively. 

3. RESULTS 

3.1. Molecular emission 

Figs. □ to a show spectral maps of the CS(2-^1), CS(3^2), C34S(3-^2), SiO(2^1), 
SiO(3^2), CH30H(3fc — > 2^), and CH30H(2fc 1^) line emission observed toward 
G324.201+0.119, G328. 307+0.423, G329. 337+0.147 and G330.949-0.174, respectively. The 
spacing is 30" for all transitions, except CS(2^1) for which the spacing is 45". Offsets 
are from the reference position given in columns 3 and 4 of Table [H Table H] gives the 
parameters of the line emission observed at the peak position (cols. 2 to 4) and of the 
spatially averaged emission (cols. 5 to 7), determined from Gaussian fits to the respective 
line profiles. 

Carbon monosulfide was detected in all observed transitions towards all sources. The 
CS(2— *1) and CS(3— >2) line profiles indicate the presence of two emission components: (i) 
a bright, spatially extended, and relatively narrow velocity feature; and (ii) a weaker and 
broader velocity feature, which is usually spatially confined to the central region. The line 
widths of the narrow component emission, which most likely arises from the dense ambient 
gas, are typically about 5 km s~^ , a value considerably larger than the thermal width 
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expected for kinetic temperatures of ~40 K, of ~ 0.2 km , implying that the hne widths 
are dominated by non-thermal supersonic motions. Whether the non-thermal motions are 
due to the effect of stellar winds or are intrinsic to the initial conditions of the formation 
process of massive stars is still an open question. 

The velocity range of the broad wing emission (full width at zero power) in different 
transitions are given in Table [51 The broad emission, which has typically full widths at zero 
power of ~ 26 km s^^ , is mainly detected toward the peak position of the massive cores. 
Although this emission most likely arises from outflowing gas, the coarse angular resolution 
(~ 30") of our observations does not allow to determine its spatial distribution and degree 
of coUimation. The nature of the high velocity gas giving rise to the pedestal feature is thus 
difficult to assess. Higher angular resolution observations are needed to investigate whether 
the high velocity emission arises from a coUimated bipolar outflow driven by the recently 
formed massive stars or from an spherically symmetric structure possible driven either by 
stellar winds or by the expansion of the HII region onto the ambient medium. 

Silicon monoxide was detected toward the four observed high-mass star forming regions, 
the strongest emission being observed toward the G330. 949-0. 174 core. SiO lines are tracers 
of both the warm ambient gas surrounding young luminous stars and of molecular outflows 
associated with newly formed massive stars. Harju et al. (1998) have already shown that 
SiO emission is relatively easy to detect toward regions of high-mass star formation, with 
the detection rate increasing with the FIR luminosity. The profiles of the SiO emission 
show a broad component of high velocity gas superposed on a narrower component, whose 
central velocity is close to the velocity of the ambient gas. The full widths at zero power 
(emission above a 3a level) of the broad component range from 17.4 to 36.6 km s~^ , with 
a median of 26 km s^^ . The relative contribution to the observed intensity is different 
from source to source. In G324.201 the SiO emission from the broad component dominates 
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the spectrum, whereas in G330.949 the narrower emission is stronger. We suggest that, 
on the one hand, the narrow SiO emission arises from the warm ambient gas surrounding 
the recently formed massive stars, and that this component resuhs from the evaporation 
of sihcon compounds from grain mantles produced by the stellar radiation. On the other 
hand, the broad SiO emission most likely indicates mass outflows from a star or stars 
embedded within the massive core, and that this component results from grain disruption 
by shock waves. Hence, the SiO emission traces both high kinetic temperature ambient gas 
and outflowing shocked gas. 

The shape of the profiles are different from molecule to molecule depending on the 
relative intensities between the narrow and broad components. To illustrate in more detail 
the differences between the line profiles. Figs. [5] to [8] show the spectra of all molecular 
lines, except methanol, observed at the peak position of each region. The CS(5^4) spectra 
corresponds to the average spectra over the 30"x30" mapped region, and therefore it is 
directly comparable to the peak spectra observed in the other transitions with ~30" beams. 

Fig. [HI shows contour maps of the velocity integrated ambient gas emission in the 
CS(2— i>l) line toward all four regions. The 50% contour level of the emission is marked with 
dark lines. In all cases the emission, at the resolution of 52", arises form a single source 
exhibiting a simple morphology. The deconvolved major and minor FWHM angular sizes 
determined by fitting a Gaussian profile to the emission are given in §3.3. 

3.2. Dust continuum 

Fig. flOl presents maps of the 1.2-mm continuum emission, within regions of typically 4' 
in size, towards the four observed sources. Emission was detected toward all four regions. 
The extent and morphology of the 1.2-mm dust continuum emission is similar to that of 
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the CS(2 —>■ 1) line emission, indicating that these two probes trace the same physical 
conditions. The observed parameters of the 1.2-mm sources are given in Table O Cols. (2) 
and (3) give their peak position. Cols. (4) and (5) give, respectively, the peak flux density 
and the total flux density, the later measured directly from the maps using the AIPS tasks 
IMEAN. Col. (6) gives the deconvolved major and minor FWHM angular sizes determined 
by fitting a single Gaussian profile to the whole spatial distribution, except for G329.337 for 
which two Gaussian components were fitted. 

In three of the four regions the morphology of the dust emission can be described 
as arising from a bright compact peak surrounded by an extended envelope of weaker 
emission, which is characteristic of sources with centrally condensed density profiles. This 
is illustrated in Fig. [TT] which shows slices of the observed 1.2-mm intensity across the 
two sources with nearly circular morphology (G324.201 and G330.949). We find that the 
observed radial intensity profiles can be well fitted with single power-law intensity profiles 
of the form / oc r^°, where r is the distance from the center, properly convolved with the 
beam of 24". We derive a indices of 1.6 and 1.9 for the G324.201 and G330.949 cores, 
respectively. For the G328.307 core an a index of 1.5 has been already derived by Garay et 
al. (2007). 



3.3. Individual sources 

Three of the high-mass star forming regions investigated here, G324.201, G329.337 
and G330.949, are in the list of ~2000 massive young stellar objects (MYSO) candidates 
compiled by the RMS survey based on a colour selection criteria of MSX point sources 
(Hoare et al. 2005). Radio continuum observations have, however, shown that all four 
Norma cores are associated with compact radio sources, indicating that they are in a more 
advanced stage of evolution in which an UCHII region has already developed (Walsh et al. 
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1998; Garay et al. 2006; Urquhart et al. 2007). In all cases the compact radio sources are 
found projected at the peak position of the 1.2-mm dust continuum emission, suggesting 
that massive stars are located at the center of the massive and dense cores. In what follows 
we discuss the characteristics of the molecular gas and dust emission, as well as of the radio 
emission, toward each of the IRAS sources, individually. 

G324.201+0.119.- The G324.201+0.119 massive core is associated with the IRAS source 
15290-5546, which has a total luminosity of 5.9 x IO^Cq. The CS profiles across the core 
show emission from a nearly Gaussian component, with central velocities of ~ —88.4 
km s~^ and FWHM widths of ~ 4.0 km s~^ . The spatial distribution of the emission in the 
CS(2 — s> 1) line from the ambient gas is slightly elongated, with FWHM deconvolved major 
and minor axis of 43" and 30", respectively. The morphology of the 1.2-mm dust emission 
also shows the presence of a single structure with FWHM deconvolved major and minor 
axis of 35" and 28", respectively. A power-law fit to the 1.2-mm intensity profile gives an 
index of 1.6 (see Fig. [TTj) . The total fiux density of the region at 1.2 mm is 17.3 Jy. 

Toward the center of the core the Gaussian component is superposed on a broader and 
weaker component. In the CS(3 — > 2) line the wing emission extends up to a fiow velocity 
of —13.1 km s^^ in the blue side and 11.4 km s^^ in the red side. The fiow velocity, Vfiow, 
is defined as Vfiow = v — Vo, where Vo is the ambient cloud velocity. The SiO line profiles 
are, on the other hand, dominated by the emission from the broad component, covering a 
fiow velocity range from —12.6 to 16.8 km s~^ . 

Radio continuum observations towards G324. 201+0. 119 show the presence of two 
compact radio sources: a cometary-like component, with an angular diameter of ~4", 
located at the center of the core and an unresolved component located ~6" north of the 
former (Urquhart et al. 2007). If excited by single stars, the rate of UV photons needed to 
produce the cometary and unresolved regions of ionized gas implies the presence of 07 and 
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08.5 ZAMS stars, respectively. 

G328.307+0.423.- The G328.307+0.423 massive core is associated with the IRAS source 
15502-5302, which has a total luminosity of 1.1 x IO^Cq. The spatial distribution of the 
CS(2 1) emission from the ambient gas shows an elongated morphology, with FWHM 
deconvolved major and minor axis of 65" and 37", respectively. The elongated morphology 
is more clearly seen in the observed spatial distribution of the 1.2-mm emission (see Fig. [TOj) . 
which has FWHM deconvolved major and minor axis of 43" and 26", respectively. The 
total flux density of the region at 1.2-mm is 24.1 Jy. 

The C0(1 0) spectra observed at the center of the core show a double peaked profile, 
with a strong blueshifted peak at the velocity of —93.9 km s^^ and a weaker redshifted 
peak at the velocity of —88.5 km s~^ . On the other hand, the spectra of the optically thin 
C^^0(1 0) line show a single Gaussian profile with a peak velocity of -92.6 km s~^ . 
These spectral characteristics indicate the presence of infall motions (Leung & Brown 1977, 
Walker, Narayanan, & Boss 1994, Myers et al. 1996). Using expression (9) of Myers et 
al. (1996) we derive an infall velocity of ~ 0.3 km s~^ . The profiles of the CS and C^^S 
emission observed across the core are generally asymmetric, showing a peak at velocities 
of ~ —92.8 km s~^ and a shoulder toward redshifted velocities. The observed shapes are 
consistent with the infall motions suggested by the C0(1 — > 0) emission but for transitions 
with smaller optical depths. In addition to infalling gas, outflowing gas is clearly detected 
toward the central region. In the CS(3^2) line, blueshifted and redshifted emissions are 
detected up to flow velocities of —8.2 km s~^ and 12.2 km s^^ , respectively. The SiO 
emission, which is detected in the velocity range from —100.0 to —82.6 km s~^ , arises 
mainly from the outflow component, but a narrower component from the ambient core gas 
is also present. 

Radio continuum observations towards G328. 307-1-0. 423 show the presence of a complex 
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region of ionized gas, consisting of a handful of low brightness, extended components and a 
bright, compact component located at the center of the core (Garay et al. 2006). The flux 
density of the later component increases with frequency whereas its angular size decreases 
with frequency. These characteristics suggest the presence of a steep gradient in the electron 
density. 

G329. 337+0. 147.— The G329. 337+0. 147 massive core is associated with the IRAS source 
15567-5236, which has a total luminosity of 1.2 x IO^Cq. The CS profiles across the core 
show emission from a bright component with a central velocity of ~ —107.7 km s~^ and a 
FWHM width of ~ 4.5 km s~^ . Strong high velocity wing emission is detected toward the 
central region in all observed species (see Fig. [7]). The blueshifted and redshifted emissions 
extend up to flow velocities of —27.6 km s^^ and 24.2 km s~^ in the CO(l^O) line, —14.8 
km s~^ and 16.8 km s~^ in the CS(3— s>2) line, and —19.2 km s^^ and 17.4 km s^^ in 
the SiO(3^2) line. 

The CS(2 — > 1) emission from the ambient gas arises from a structure that is slightly 
elongated in the SE direction, with deconvolved FWHM major and minor axis of 44" and 
39", respectively. The peak of the CS emission is coincident with a methanol maser (Walsh 
et al. 1997). The spatial distribution of the silicon monoxide and methanol emission 
show notable differences with respect to that of the emission in carbon monosulfide. The 
morphology of the SiO and CH3OH emission shows a double peaked structure, with the 
main peak being displaced from the peak of the CS emission by ~ 60" toward the southeast 
(see Fig. [T^ . The secondary SiO peak coincides with the peak of the CS emission. The 
double peaked morphology is also clearly seen in the map of the dust continuum emission. 

Radio continuum observations towards G329. 337+0. 147 show the presence of a single 
source, with an angular diameter of 6.3 x 4.8", exhibiting an irregular multi-peaked 
morphology (Walsh et al. 1998; Urquhart et al. 2007). If excited by a single star, the rate 
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of UV photons required to ionize this region imphes the presence of an 06.5 ZAMS star. 

G330 .949-0. 174.- The G330.949-0.174 massive core is associated with the IRAS source 
16060-5146, which has a total luminosity of 1.0 x IQ^Cq. The CS(2— i>l) line exhibit an 
asymmetric profile, with a peak at the velocity of —92.1 km s~^ and a shoulder toward 
redshifted velocities. The CS(3^2) and CS(5— s>4) profiles are double peaked, with peaks at 
~ —93.0 km s~^ and ~ —89.1 km s~^ . The blue peak is stronger than the red peak by 
a factor of ~ 1.5. The profiles of the C^^S and SiO emission are, on the other hand, more 
symmetric, with peak velocities lying between the two velocity peaks of the CS emission. 
This result indicates that the complex CS profiles are produced by optical depth effects 
rather than being due to two velocity components along the line of sight. 

The presence of a high velocity gas toward the central region is clearly indicated by 
the profiles of the CS and SiO emission. In CO lines the outflow emission is blended 
with emission from less dense molecular clouds in the line of sight. The emission in the 
CS(3 — » 2) line extends up to a velocity of —75.7 km s^^ in the red side and up to a 
velocity of —106.0 km s~^ in the blue side. 

The spatial distribution of the CS(2 1) emission from the ambient gas is nearly 
circular, with FWHM deconvolved major and minor axis of 43". The circular morphology 
is also seen in the 1.2-mm dust emission, which exhibits a single structure with FWHM 
deconvolved major and minor axis of 33". The total flux density of the region at 1.2 mm 
is 47.2 Jy. A power-law fit to the 1.2-mm intensity profile gives an index of 1.9 ± 0.1 (see 
Fig. [TT]) . The peak of the molecular and dust emissions are displaced about 30" east of the 
IRAS position, but are coincident with the position of a methanol maser and a compact 
radio continuum source (Walsh et al. 1997, 1998). The radio continuum observations of 
Urquhart et al. (2007) towards G330. 949-0. 174 show the presence of two compact sources: 
a bright cometary-like component, with an angular diameter of ~3", located near the 
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center of the core and an unresolved, weak component, located ~12" northeast of the 
former. If excited by single stars, the rate of UV photons needed to produce the cometary 
and unresolved regions of ionized gas implies the presence of 07 and BO ZAMS stars, 
respectively. 
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4. DISCUSSION 

4.1. Spectral Energy Distributions 

Fig. [13] shows, for each of the 4 sources, the spectral energy distribution (SED) from 
12 fim to 1.2 mm. In this wavelength range the emission is mainly due to thermal dust 
emission. The flux densities at 12, 25, 60 and 100 fim correspond to IRAS data, and those 
at 8.3, 12.1, 14.7 and 21.3 fim correspond to Midcourse Space Experiment (MSX) data 
(Price et al. 2001). The IRAS and MSX fluxes were obtained from the IPAC database. 

We fitted the SED with modified blackbody functions of the form ByiT^) [1 — exp(— Tj,)] Vl 
where is the dust optical depth, B^(Td) is the Planck function at the dust temperature 
Trf, and is the solid angle subtended by the dust emitting region. The opacity was 
assumed to vary with frequency as u^, i.e. t^, = {u/uo)^, where Uo is the frequency at which 
the optical depth is unity. A single temperature model produced poor fits, underestimating 
the emission observed at wavelengths smaller than 25/im, and therefore we used a model 
with two temperature components. Since the regions are likely to present temperature 
gradients, this is indeed a coarse simplification. However, this model allows us to determine 
the average dust temperatures representative of each wavelength range. The parameters 
derived from the fits, namely dust temperature, opacity power-law index, wavelength at 
which the opacity is unity, and angular size (assuming a Gaussian flux distribution), for the 
colder dust component are given in Fig. [131 



4.2. Parameters derived from dust observations 

The parameters derived from the 1.2-mm observations are summarized in Table [71 The 
dust temperatures given in col. (2) correspond to the colder temperature determined from 
fits to the spectral energy distribution (SED). The radius given in col. (3) were computed 
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from the geometric mean of the deconvolved major and minor angular sizes obtained from 
single Gaussian fits to the observed spatial structure. The masses of the cores, given in 
col. (4), were derived using the expression, (e.g., Chini, Kriigel, & Wargau 1987) 

Rdg 1 . 2mm Bi . 2mm {Td) 

where Si_2mm is the observed flux density at 1.2mm, Ki_2mm is the mass absorption coefficient 
of dust with a value of 1 cm^ g~^ (Ossenkopf & Henning 1994), R^g is the dust-to-gas mass 
ratio (assuming 10% He) with a value of 0.01, and B{Td) is the Planck function at the 
dust temperature T^. This expression implicitly assumes that the observed emission at 1.2 
mm corresponds to optically thin thermal dust emission and that the source is isothermal. 
The masses derived in this way range from 5 x 10^ to 1 x 10^ Mq. Cols. (5) and (6) give, 
respectively, the average molecular densities and average column densities derived from the 
masses and radius assuming that the cores have uniform densities. Clearly, this is a rough 
simplification, since as discussed in §3.2 the massive and dense cores are likely to have steep 
density gradients. Finally, the continuum optical depth at 1.2- mm is given in col. (7). 



4.3. Parameters derived from line observations 

J,.. 3.1. Optical depths and column densities of CS and CO 

The optical depth in the CS(3— i^2) line, r^i' , can be derived from the ratio of observed 
main beam brightness temperatures in the CS(3^2) and C^^S(3^2) lines using the 
expression, 

1 - exp(-r3^//r) ^ [Jg'{T^.) - Jg'{T,g)] T^,{C^'S) 
1 - exp(-r3^/) [Jg^'^iT^^) - J3S"^(n,)] T^,{CS) ' ^ ^ 

where the subscript 32 refers to the 3^2 transition, T^x is the excitation temperature of 

the transition, T^g is the background temperature, 

hv 1 



k ex.p{hiy/kT) - 1 
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and r is the /t§2 ^ optical depth ratio, given by 

I p f exp(-12/ii?^7fcTe..) ) f e^h^^' /kT,^) - 1 1 {kTjhB^''^ + 1/3) 
I /x'^^^^^ P 1 exp(-12/ii?t^^^^^/A;Te,.) J 1 eMhi^^"' /kTe.) " 1 J {kTeJhB^^ + 1/3) ' 

(3) 

where a is the [C'^^S/C'^^S] isotopic abundance ratio, and /x and B are the permanent 
dipole moment and the rotational constant of the molecule, respectively. If the excitation 
temperature is known, the above expressions allow the optical depth to be determined. We 
avoid the usual approximation that ^ 1, and solve equation ([2]) using an interpolation 
procedure, assuming a [C^^S/C^^S] abundance ratio of 22.5 (Blake et al. 1994). 

The optical depth in the C0(1 — > 0) line, be derived from the ratio of 

observed main beam brightness temperatures in the C0(1 —* 0) and ^^C0(1 — > 0) lines in a 
similar manner as described above (see Bourke et al. 1997, for the appropriate expressions). 
We assumed a [CO/^'^CO] abundance ratio of 55, corresponding to the average value across 
the Galactic disk (Wannier et al. 1982). Further, the observations of the J=l line 
emission in the three isotopic species of carbon monoxide (CO, ^^CO and C^^O) allow to 
determine the [^^COj/fC^'^O] isotopic abundance ratio. Using a [CO]/[^'^CO] abundance 
ratio of 55, we derive [^^CO]/[C^^O] abundance ratios in the range 8 to 10. 

The peak optical depths at the peak position of the massive cores in the CS(3— s>2), 
C34S(3^2), CO(l^O) and 1^00(1^0) lines, are given in cols. 2 to 5 of Tabled They were 
computed assuming Tex = 40 K. We find that the emission is optically thick in the CS lines 
and optically thin in the C^^S line for all cores. We note that the dependence of the optical 
depth with excitation temperature is weak. For instance, assuming T^.^ = 100 K the derived 
optical depths are only 0.4% larger. 

The total column density of a linear, rigid rotor, molecule can be derived from the 
optical depth and excitation temperature of a rotational transition at frequency z/, using 
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the expression (e.g. Garden et al. 1991) 

SttV {J + l)hB 1 -exp(-/ii^/A;Te,) y ^ ^ ' 
where J is the rotational quantum number of the lower state. This expression assumes that 
all energy levels are populated according to local thermodynamic equilibrium (LTE) at the 
temperature Tg^- In particular the total column density of the CS molecule is given, in 
terms of the opacity and excitation temperature of the 3^2 transition, by 

iV(CS) = 5.93 X 10^^ (Te. + 0.39)-^^^^^^^^^^^ [ rg' dv , (5) 

1 - exp(-7.05/Te^) J 

where v is measured in km s^^ . Col. 2 of Table [9] gives the derived CS column density 
through the center of the cores. The CO column density through the center of the cores 
are given in col. 3 of Table [91 They were determined in the same way as the CS column 
densities using the [CO,^^CO] pair of observations. 

From the column density in CS, N(CS), we can estimate the mass of the cores using 
the expression, 



M = N{CS) 



Ho 



fi^ niH TcR^ , (6) 



_CS_ 

where [^] is the H2 to CS abundance ratio, jjm is the mean molecular weight, and R is the 
core radius. The radii, given in column 2 of Table [TOl were computed from the geometric 
mean of the deconvolved CS(2— 5>1) angular sizes using the distances given in Tabled! The 
core masses computed from expression (6), using [^] = 3 x 10'' (see §4.3), are given in 
col. 3 of Table [TOl They range from 4.3 x 10'^ to 1.0 x 10^ Mq. These masses are in good 
agreement with those derived from the dust continuum observations. 



4-3.2. Rotational temperatures and column densities of methanol and silicon monoxide 



The rotational temperature. Trot, and the total column density, Nt, of methanol and 
silicon monoxide can be computed using a rotational diagram analysis (see e.g., Linke, 
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Frerking & Thaddeus 1978; Blake et al. 1987), with the assumptions of optically thin 
conditions and local thermodynamical equilibrium (LTE), which relates the integrated line 
intensity, rotational temperature, and column density via 

3A; / T^hdv ^ ^ Nt ( (7) 

87t''^IJ,'^uS Qu Q{Trot) \ kTrot ) 

where /x, and 5* are the transition dipole moment, frequency, and line strength of the 
transition, respectively, j Tmbdv is the velocity integrated main beam brightness, obtained 
directly from the observations, Eu is the upper state energy, and Q{Trot) is the rotational 
partition function. 

Fig. [H] shows rotational diagrams for the CH3OH emission observed at the peak 
position of the massive cores. The column densities and rotational temperatures, derived 
from a linear least squares fit to the data, are given in cols. 4 and 6 of Table [91 The 
rotational temperatures of methanol range from 12.0 K to 22.6 K. These temperatures are 
significantly smaller than the temperatures derived from the dust observations, of ~ 40K, 
indicating that the methanol populations are sub-thermally excited (Bachiller et al. 1995; 
Avery & Chiao 1996). 

From rotational diagram analysis of the SiO emission observed at the peak position 
of the massive cores we derived column densities and rotational temperatures of both the 
ambient gas, given in cols. 5 and 7 of Table [9l and flowing molecular gas, given in cols. 8 
and 9. The rotational temperatures of the flowing gas are typically two times larger than 
those of the ambient gas. This is illustrated in Fig. [T51 which presents a rotational diagram 
for the SiO emission observed at the peak position of the G330.949 massive core. In this 
case the gas giving rise to the narrow emission has a rotational temperature of 11.2±0.6 K 
and a column density of 4.1 x 10^^ cm^^, whereas the outflowing gas producing the broad 
emission has a rotational temperature of 24.9±2.4 K and a column density of 2.6 x 10^^ 

_2 

cm . 
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4-3.3. Molecular abundances 

The chemical composition of a gas cloud is usually characterized by the fractional 
abundance of molecules relative to molecular hydrogen, the main constituent of the 
interstellar gas. The abundance of H2 is however not directly measured. Our observations 
allow us, however, to directly determine abundances relative to CO. We compute the 
abundance of species X relative to CO, [X/CO] (where X = CS, SiO, or CH3OH), as the 
ratio of the column density of species X and the CO column density in the same velocity 
range. For SiO and CH3OH we use the column densities derived from rotational analysis 
(see §x.x) whereas for CS and CO we use the column densities derived from the observations 
of their main and isotopic species (see §4.3.1). Table [TT] summarizes the abundances of CS, 
SiO and CH3OH relative to CO. To convert [X/CO] to fractional abundances relative to H2 
we will assume [CO/H2]=10~^. 

Silicon monoxide. We find that the [SiO/CO] abundance ratio in the ambient gas range 
from 3 X lO""^ to 1 x 10~^. Assuming that [CO/H2] is 10"'', then the [SiO/H2] abundance 
ratio in the ambient gas ranges between 3 x 10^^^ and 1 x 10^^°. Chemical models predict 
that both high temperature and shocks can produce enhancements in the gas phase SiO 
abundance. The abundance of SiO in the ambient gas is likely due to the gas-phase 
high temperature chemistry that takes place in the hot regions surrounding the recently 
formed massive stars. Models of gas-phase chemistry of hot molecular cores that takes into 
account grain mantle evaporation by stellar and/or shock radiation (MacKay 1995), predict 
[SiO/H2] abundances of 3 x 10"^^ for gas kinetic temperatures of 100 K and densities of 
2 X 10^ cm-^ 

We also find that the [SiO/CO] abundance ratios in the outflowing gas are similar to 
those derived in the ambient gas. The abundance of SiO in the high velocity gas is likely to 
be due to gas phase chemistry behind a shock (Iglesias & Silk 1978; Hartquist et al. 1980; 
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Mitchell & Deveau 1983). The enhancement of the SiO abundance in the high velocity 
material is probably the result of grain destruction by shocks which releases silicon into the 
gas phase, allowing SiO to form (Scab & ShuU 1983; Schilke et al. 1997; Caselli et al. 1997). 

Methanol. We find that the [CH3OH/CO] abundance ratio in the ambient gas is similar 
for all observed cores, having a mean value of 6 x 10~^. Assuming an [H2/CO] of 10'', it 
imphes that the [CH3OH/H2] abundance ratio in the ambient gas is typically 6 x 10~^. At 
temperatures smaller than 100 K, chemistry models predict that the production of CH3OH 
in the gas phase yield abundances of only 10~^^ relative to H2 (Lcc ct al. 1996). The 
high abundance of methanol in the cores suggests that the radiation from their embedded 
luminous stars are producing a substantial evaporation of icy grain mantles. Van der Tak et 
al. (2000a) found a [CH3OH/H2] ratio of a few 10~^ for hot cores and ~ 10~^ for the cores 
with the lower temperatures. The values obtained for the cores investigated here are in 
between the values found for dark clouds and hot cores, and may correspond to an average 
of the values appropriate for the cold envelope mass and hot core mass. 

Carbon monosulfide. We find that the [CS/CO] abundance ratio in the ambient gas is 
similar for all observed cores, having a mean value of 4 x 10"^. Assuming an [H2/CO] of 
lO"', it implies that the [H2/CS] abundance ratio in the ambient gas is typically 3 x 10^. 
This value is typical of molecular cores associated with UC HII regions and hot cores (van 
der Tak et al. 2000b). 

4.3.4- Structure of the cores 

In the previous analysis we have implicitly assumed that the density and temperature 
arc uniform within the massive cores. This is obviously a simplification, since the physical 
conditions are likely to vary with distance from the center of the core. In fact, the 
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intensity profiles of the dust continuum emission indicate that the cores are highly centrally 
condensed. 

Assuming that cores have density and temperature radial distributions following 
power laws, then for optically thin dust emission the intensity index a is related to the 
density index p (n oc r~^) and the temperature index q [T (x r~^), by the expression 
a — p + Qq — 1 + €f (Adams 1991; Motte & Andre 2001), where Q is a temperature and 
frequency correction factor with a value of ~ 1.2 at 1.2-mm and 30 K (Beuther et al. 2002) 
and 6/ is a correction factor to take into account the finite size of the cores. The massive 
and dense cores investigated here are heated by a cluster of luminous sources located at 
their central position (Chavarria et al. 2009; Paper II). For this type of objects van der Tak 
et al. (2000b) found that the temperature decreases with distance following a power-law 
with an index of 0.4. Adopting this value for q and a value of 0.1 for e/, we then infer, 
using the above expression, that massive and dense cores have density distributions with 
power-law indices in the range 1.9 — 2.3. Several dust continuum studies have already 
shown the presence of density gradients within massive and dense cores (e.g., van der Tak 
et al. 2000b, Mueller et al. 2002, Beuther et al. 2002, WiUiams et al. 2005). In particular, 
Beuther et al. (2002) observed a large sample of high-mass protostellar candidates and 
found that the steeper indices are associated with the more luminous and more massive 
objects. The steep indices derived for the Norma cores are characteristics of the subgroup 
of strong molecular sources within Beuther's sample, which have a mean density index of 
1.9. The steep profiles are thought to indicate objects in the collapse and accretion phase. 

4.4. Segregation within the massive and dense cores 

Given the large mass of the massive and dense cores investigated here, of typically 
5 X 10"^ Mq, their overall collapse process is likely to produce a protostar cluster. The 
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questions of how an individual massive star forms and how the cluster forms are therefore 
closely related. The derived properties of massive and dense cores are then important 
boundary conditions to be taken into account for models and simulations. 

For dynamical considerations, the relevant mass of the cores is their total mass, 
including gas and stars. We have already shown that the mass of gas in the cores 
determined from two independent methods, one based on dust observations and the other 
on molecular line observations, are in very good agreement. The determination of the total 
mass of the cores is not straightforward and is usually estimated assuming that the cores are 
in virial equilibrium. Masses computed under this assumption, using the observed velocity 
dispersion and size of the CS(2— emission, are given in column 4 of Table [TOl The virial 
masses are similar, within the errors, to the gas masses, implying that the bulk of the core 
mass is in the form of molecular gas. We note that the virial masses are in average a factor 
of 1.5 smaller than the dust or molecular masses which might also indicate that the cores 
do not have enough support and may be contracting. 

In paper II we show that the massive cores investigated here are associated with 
clusters of embedded stars which exhibit clear mass segregation. The high-mass stars are 
found concentrated near the center of the core while intermediate mass stars are spread 
over much larger distances from the core center. Garay et al. (2006) have already reported 
that massive stars are typically located near the center of the massive and dense cores. The 
question arises as to whether the massive stars are formed in situ or if their spatial location 
is a result of dynamical effects. 

The characteristics of the cores derived here, namely their high densities, steep density 
profiles and sizes, together with the relatively low fraction of the total mass in the form 
of stars, suggest that dynamical friction (also called gravitational drag) produced by the 
gaseous background onto the stars may play an important role in the orbital evolution of 
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stars (Ostriker 1999, Sanchez-Salcedo & Brandenburg 2001, Kim & Kim 2007). The role of 
gravitational drag on a gaseous background in the evolution and migration of stars within 
dense gaseous star forming clouds has been investigated using numerical smooth particle 
hydrodynamics simulations (Escala et al. 2003, 2004). Simulations have been carried out 
with initial conditions similar to those of the dense cores investigated here, namely steep 
cloud density profile (n(r) oc r^^-^) and perturber having masses of 1 % of the total gas mass 
(Mstar ~ 20 Mq). For a typical cloud mass of 4 x 10^ Mq and initial distance of the stars to 
the center of the cloud of 0.4 pc, the migration timescale for high-mass stars is 4.5 x 10^ yr. 
This value is smaller than the estimated ages of the young clusters within the cores of ~ 
1 Myr (see paper II). Since the migration timescales for dynamical friction are proportional 
to the inverse of the perturber's mass (Binney & Tremaine 1987), intermediate mass stars 
have migration timescales typically a factor ~ 10 longer than those of high mass stars and 
greater than the estimated age of the clusters. Therefore, gravitational drag is considerably 
less effective for intermediate mass stars than for high mass stars. This suggest that in the 
young embedded stellar clusters only the massive stars have been able to efficiently migrate 
towards the center due to gravitational drag. We conclude that dynamical interactions of 
cluster members with the gaseous component of massive and dense cores can explain the 
observed mass segregation, indicating that the gas play an important role in the dynamic 
of young cluster members and originating the observed mass segregation. 

We can not rule out, however, the hypothesis that the cradles of massive stars are 
formed by direct accretion at the center of centrally condensed massive and dense cores 
(McKee & Tan 2003). In their turbulent and pressurized dense core accretion model, the 
collapse of a massive and dense core is likely to produce the birth of a stellar cluster, with 
most of the mass going into relatively low-mass stars. The high-mass stars are formed 
preferentially at the center of the core, where the pressure is the highest, and in short time 
scales of - 10^ yrs (Osorio et al. 1999; McKee & Tan 2002). 



5. SUMMARY 



We undertook molecular line observations in several species and transitions and dust 
continuum observations, both made using SEST, towards four young high-mass star forming 
regions associated with highly luminous IRAS point sources with colors of UC HII regions 
and CS(2— >1) emission. These are thought to be massive star forming regions in early stages 
of evolution. The objectives were to determine the characteristics and physical properties 
of the dust clouds in which high-mass stars form. Our main results and conclusions are 
summarized as follows. 

We find that the luminous star forming regions in Norma are associated with molecular 
gas and dust cores with radii of typically 0.5 pc, masses of ~ 5 x 10^ Mq, molecular 
hydrogen densities of typically ~ 2 x 10^ cm~^, column densities of ~ 5 x 10^^ cm~^, and 
dust temperatures of ~ 40 K. All the derived physical parameters of the Norma cores, 
except the last, are similar to those of cores harboring UC HII regions (Hunter et al. 2000; 
Faiindez et al. 2004). For a sample of about 150 IRAS point sources with colours of UC 
HII regions and average luminosity of 2.3 x 10^ L©, Faundez et al. (2004) derived an 
average mass of 5 x 10'^ Mq, an average size of 0.4 pc, and an average density of 2 x 10^ 
cm~^. The dust temperature of the Norma cores are however larger than of the Faundez et 
al. cores, for which they find an average temperature of 32 K. The larger values derived for 
the Norma cores are consistent with them being centrally heated by more luminous objects 
(average luminosity of 9.7 x 10^ Lq). 

We find that the observed radial intensity profiles of the dust continuum emission from 
the Norma cores are well fitted with single power-law profiles, with intensity indices in the 
range 1.5 to 1.9. This in turn indicates that the cores are highly centrally condensed, having 
radial density profiles with power-law density indices in the range 1.9 — 2.3. These steep 
profiles are thought to indicate objects in the collapse and accretion phase of evolution. 
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We find that under the physical conditions of the Norma cores, high-mass stars can 
migrate towards the central region due to gravitational drag in time scales of 5 x 10^ yr, 
smaller than the estimated ages of the clusters, providing an explanation for the observed 
stellar mass segregation within the cores (Paper II). 

The authors gratefully acknowledge support from CONICYT through projects 
FONDAP No. 15010003 and BASAL PFB-06. 
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Table 1. OBSERVED SOURCES 



Source 


IRAS 


a (2000) 


5(2000) 


D 












(kpc) 


( L .) 


G324.201+0.119 


15290-5546 


15^32"^53!2 


-55°56'13" 


6.0 


5.9 X 10^ 


G328.307+0.423 


15502-5302 


15 54 06.3 


-53 11 38 


5.6 


1.1 X 10^ 


G329.337+0.147 


15567-5236 


16 00 33.4 


-52 44 45 


7.0 


1.2 X 10^ 


G330.949-0.174 


16060-5146 


16 09 49.1 


-51 54 45 


5.4 


1.0 X 10^ 



Tabic 2. MOLECULAR LINES: OBSERVATIONAL PARAMETERS 

Molecule Transition Frequency Beam Size 77 Ai> 

(MHz) (FWHM ") ( km s"^ ) 



CS 



Q34S 

SiO 

CH30H 

CO 
i^CO 

C^SQ 



J = 2 
J = 3 
J- 5 
J = 3 
J = 2 
J = 3 
</fc = 2o — 

Jk — 3-1 - 

J = 1 
J = 1 
J = 1 



1 

2 
4 
2 
1 
2 

lo A+ 
> 2_i E 









97980.968 


52 


0.73 


0.132 


146969.049 


34 


0.66 


0.085 


244935.606 


22 


0.45 


0.053 


144617.147 


34 


0.66 


0.087 


86846.998 


57 


0.75 


0.144 


130268.702 


40 


0.68 


0.098 


96741.42 


52 


0.73 


0.129 


145097.470 


34 


0.66 


0.086 


115271.204 


45 


0.70 


0.111 


110201.353 


47 


0.71 


0.116 


109782.160 


47 


0.71 


0.116 
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Table 3. INTEGRATION TIMES AND SENSITIVITIES 



Line G324.201 G328.307 G329.337 G330.949 



CS(2 ^ 1) 
CS(3 ^ 2) 
CS(5 ^ 4) 
C^S{3 2) 
SiO(2 ^ 1) 
SiO(3 ^ 2) 
CH30H(2fe ^ U) 
CH30H(3fe ^ 2k) 
C0(1 ^0) 
i3C0(l ^ 0) 
Ci«0(l ^ 0) 



tin (T tjyj (T tjjj (T 

(m) (K) (m) (K) (m) (K) (m) (K) 



3 


0.12 


3 


0.14 


6 


0.064 


3 


0.094 


4 


0.065 


4 


0.071 


4 


0.065 


4 


0.057 


5 


0.47 


3 


0.11 


6 


0.41 


5 


0.44 


4 


0.074 


4 


0.068 


4 


0.068 


4 


0.049 


20 


0.018 


12 


0.028 


12 


0.027 


8 


0.026 


12 


0.033 


22 


0.016 


22 


0.021 


4 


0.035 


8 


0.037 






4 


0.054 


4 


0.038 


8 


0.043 


8 


0.045 


8 


0.044 


4 


0.049 


1 


0.276 


1 


0.276 


1 


0.254 


1 


0.254 


3 


0.072 


3 


0.071 


3 


0.068 


3 


0.076 


6 


0.057 


6 


0.055 


5 


0.060 


5 


0.063 



^1(7 rms noise in antenna temperature. 
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Table 4. MOLECULAR LINES: OBSERVED PARAMETERS 



Line 






Peak ^ 






Average ^ 










V 


Av 




V 


Av 






(K) 


( km ) 


( km ) 


(K) 


( km ) 


( km ) 




G324.201+0.119 


CS(2^ 1) 




1.79±0.12 


-88.44±.04 


4.1±0.3 


0.537±.02 


-88.16 ± .03 


4.6 ±0.1 


CS(3 ^ 2) 




2.17±0.07 


-88.52±.01 


4.1±0.1 


1.31±.02 


-88.34 ± .09 


4.4 ±0.1 


CS(5 ^ 4) 




2.98±0.47 


-87.87±.05 


4.4±0.1 


2.33±0.16 


-88.33 ± .04 


3.8 ±0.1 


0=^^8(3 ^ 2) 




0.745±0.074 


-88.59±.05 


3.6±0.1 


0.276 ± .02 


-88.33 ± .05 


3.4 ±0.2 


SiO(2 ^ 1) 




0.052±0.018 


-87.4±0.3 


6.1±1.1 


0.025 ± .006 


-88.74 ± .18 


3.4 ±0.5 


SiO(3 ^ 2) 










0.017 ±.011 


-88.6 ±0.3 


2.6 ±0.7 


CH30H(2o ^ 


lo)A+ 


0.380±0.037 


-88.26±.05 


4.6±0.2 


0.270 ± .012 


-88.19 ± .02 


4.2 ±0.1 


CH30H(3o ^ 


2o)A+ 


0.510±0.043 


-88.36±.03 


5.5±0.1 


0.277 ± .014 


-88.24 ± .02 


4.4 ±0.1 


G328. 307+0. 423 


CS(2 ^ 1) 






_c 




0.858 ± 0.028 


-92.65 ± .03 


6.3 ±0.1 


CS(3 ^ 2) 




2.93 ± 0.07 


-92.90±.02 


4.8±0.1 


1.84±0.02 


-92.92 ± .01 


4.3 ±0.1 


CS(5 ^ 4) 




l.lSibO.ll 


-92.82±.04 


3.6±0.2 








C34S(3 ^ 2) 




0.45 ± 0.068 


-92.49±.09 


4.4±0.2 


0.386±0.02 


-92.32 ± .03 


5.3 ±0.1 


SiO(2 ^ 1) 




0.086±0.028 


-92.24±.13 


6.5±0.4 








SiO(3 ^ 2) 




0.099 ± 0.016 


-92.75±.10 


6.1±0.3 
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Table 4 — Continued 



Line 






Peak ^ 






Average ^ 










V 


Av 




V 


Av 






(K) 


( km ) 


( km ) 


(K) 


( km ) 


( km ] 


CH30H(3o ^ 


2o)A+ 


0.37±0.045 


-91.84±.05 


5.4±0.1 


0.273 ± 0.011 


-92.71 ± .02 


4.9 ±0.1 


G329.337+0.147 


CS(2 ^ 1) 




3.34±0.06 


-107.69±.02 


4.3±0.1 


1.70±0.01 


-107.58 ±.01 


5.2 ±0.1 


CS(3 ^ 2) 




4.92±0.07 


-107.71±.01 


4.7±0.1 


1.81±0.02 


-107.57 ±.01 


4.5 ±0.1 


CS(5 ^ 4) 




2.23±0.41 


-107.73±.08 


4.5±0.3 


1.91±0.14 


-107.76 ±.04 


4.6 ±0.2 


C34S(3 ^ 2) 




1.15±0.07 


-107.72±.03 


3.6±0.1 


0.246±0.02 


-107.66 ± .02 


3.0 ±0.1 


SiO(2 ^ 1) 




0.104±0.027 


-107.67±.10 


4.2±0.2 


0.067±0.007 


-107.58 ±.05 


4.5 ±0.2 


SiO(3 ^ 2) 




0.125 ± 0.021 


-107.47±.13 


4.2±0.3 








CH30H(2o ^ 


lo)A+ 


0.502±0.054 


-107.51±.04 


3.8±0.1 


0.355±0.014 


-107.50 ± .02 


4.0 ±0.1 


CH30H(3o ^ 


2o)A+ 


0.671±0.044 


-107.48±.02 


4.8±0.1 


0.322 ± 0.012 


-107.48 ± .02 


4.5 ±0.1 


G330.949-0.174 


CS(2 ^ 1) 






_c 




0.585 ± 0.019 


-91. 16 ±.02 


7.2 ±0.1 


CS(3 ^ 2) 






_c 




1.80 ±0.016 


-91.65 ± .02 


7.2 ±0.1 


CS(5 ^ 4) 






_c 




3.05 ±0.15 


-92.12 ± .03 


9.0 ±0.1 


C3^S(3 ^ 2) 




0.624zb0.049 


-91.23±.06 


4.8±0.5 


0.270 ± 0.014 


-91.01 ± .02 


4.3 ±0.2 


SiO(2 ^ 1) 




0.501±0.026 


-90.24±.03 


6.2±0.2 


0.250 ± 0.08 


-90.3 ±0.1 


6.0 ±0.1 
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Table 4 — Continued 



Line Peak ^ Average ^ 

t;;^ V Av V Av 

(K) ( km s"^ ) ( km s'^ ) (K) ( km 8"^ ) ( km s"^ ) 

SiO(3^2) 0.461 ±0.035 -90.56±.01 6.1±0.2 _ _ _ 

CH30H(2o ^ lo)A+ 0.985 ± 0.038 -90.35±.03 6.1±0.1 0.538 ± 0.011 -90.3 ± 0.1 5.7 ± 0.1 

CH30H(3o ^ 2o)A+ 1.24 ± 0.049 -90.95±.02 7.9±0.1 0.549 ± 0.014 -90.7 ±0.1 6.9 ±0.1 

^Fits to the spectra at peak position. 
^Fits to the spatially averaged emission. 
^Non-Gaussian profile. 
'^Only broad emission detected in this line. 
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Table 5. WING EMISSION: VELOCITY EXTENT 

Source Line ^min ^max Av 

( km ) ( km ) ( km ) 

G324.201+0.119... CS(3 ^ 2) -101.5 -77.0 24.5 

SiO(2 ^ 1) -101.0 -71.6 29.4 

SiO(3^ 2) -101.4 -72.0 29.4 

G328.307+0.423... CS(3 ^ 2) -101.0 -80.6 20.4 

SiO(2 ^ 1) -100.0 -82.6 17.4 

G329.337+0.147... CS(3 ^ 2) -122.5 -90.9 31.6 

SiO(2 1) -114.0 -92.2 21.8 

SiO(3-> 2) -126.9 -90.3 36.6 

G330.949-0.174 CS(3 ^ 2) -106.0 -75.7 30.3 

SiO(2 ^ 1) -107.0 -70.4 36.6 

SiO(3^ 2) -103.4 -75.2 28.2 
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Table 6. 1.2-MM CONTINUUM EMISSION: OBSERVED PARAMETERS 



SIMBA Peak position Flux density*^ Angular size ^ 

source q;(2000) 5(2000) Peak Total (") 



(Jy/beam) (Jy) 



G324.201+0.119 


15 


32 


52.91 


-55 


56 


10.8 


5.28 


17.3 


35 


X 


28 


G328.307+0.423 


15 


54 


06.33 


-53 


11 


43.1 


6.35 


24.1 


43 


X 


26 


G329.337+0.147A 


16 


00 


33.10 


-52 


44 


46.6 


4.35 


11.7 


37 


X 


27 


G329.337+0.147B 


16 


GO 


38.68 


-52 


45 


31.7 


1.45 


4.58 


53 


X 


38 


G330.949-0.174 


16 


09 


52.41 


-51 


54 


58.5 


21.75 


47.2 


34 


X 


33 



^Errors in the flux density are dominated by the uncertainties in the flux calibration, 
of ~20%. 

^Errors in the angular sizes are typically 10%. 

Table 7. DUST EMISSION: DERIVED PARAMETERS 



SIMBA Td Radius Mass n(H2) N(H2) Ti.2mm 



source 


(K) 


(pc) 


(Mo) 


(cm' 




(cm 






G324.201+0.119 


37. 


0.45 


5.0 


X 10^ 


2.2 


X 


10^ 


4.1 


X 


1023 


0.016 


G328.307+0.423 


42. 


0.46 


5.2 


X 10^ 


2.2 


X 


10^ 


4.2 


X 


1023 


0.017 


G329.337+0.147A 


40. 


0.53 


4.1 


X 10^ 


1.1 


X 


10^ 


2.5 


X 


1023 


0.010 


G329.337+0.147B 


40. 


0.76 


1.6 


X 10^ 


1.5 


X 


10^ 


4.8 


X 


1022 


0.002 


G330.949-0.174 


37. 


0.44 


1.1 


X 10^ 


5.3 


X 


10^ 


9.6 


X 


1023 


0.038 
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Table 8. MOLECULAR LINES: OPACITIES 



Source 




Peak optical depth 






CS(3^2) 


C-^^S(3^2) 


CO(l^O) 


"CO(l^O) 


G324.201+0.119 


7.0 


0.30 


42 


0.7 


G328.307+0.423 


3.4 


0.15 


> 26 


> 0.4 


G329.337+0.147 


5.5 


0.24 


31 


0.5 


G330.949-0.174 


7.1 


0.31 


69 


1.1 



Table 9. COLUMN DENSITIES AND ROTATIONAL TEMPERATURES' 



Source 








g;is 






Oulf 


]( )\\" 


N(CS) 


N(CO) 


N(CH30H) 


N(SiO) 


T(CH30H) 


T(SiO) 


N(SiO) 


T(S 


G324.201+0.119 


7.6 X 10^^ 


1.7 X 10^9 


9.7 X 10^^ 




12.0 




1.3 X 10^3 


31 


G328.307+0.423 


5.2 X 10^^ 


1.4 X 10^^ 


9.4 X 10^"^ 


4.7 X 10^2 


10.3 


6.3 


3.8 X 10^2 


24 


G329.337+0.147 


5.8 X 10^^ 


1.4 X 10^9 


9.0 X 10^^ 


6.3 X 10^2 


17.0 


16.7 


1.8 X 10^2 


45 


G330.949-0.174 


1.6 X 10^6 


3.9 X 10^9 


2.1 X 10^5 


4.1 X 10^3 


22.6 


11.2 


2.6 X 10^3 


24 



^N: column density in cm ^; T: rotational temperature in K 
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Table 10. MOLECULAR LINES: DERIVED PARAMETERS OF CORES 



Source Radius Mass' 



(pc) ( Mo) ( Me) 



G324.201+0.119 


0.52 


4.3 


X 


10^ 


2.3 


X 


10^ 


G328.307+0.423 


0.66 


4.6 


X 


10=^ 


5.5 


X 


10^ 


G329.337+0.147 


0.71 


6.0 


X 


10-^ 


4.0 


X 


10^ 


G330. 949-0. 174 


0.56 


1.0 


X 


101 


6.1 


X 


10'^ 



^Derived from CS(2— >1) observations. 



Table 11. ABUNDANCE RATIOS 



Source CS/CO SiO/CO CH3OH/CO 



G324.201+0.119 4.5 x 10"^ 

G328.307+0.423 3.7 x 10"^ 

G329.337+0.147 4.1 x 10"^ 

G330.949-0.174 4.1 x 10"^ 



— 5.7 X 10-^ 

3.3x10-^ 6.7x10^5 

4.5 X 10"^ 6.4 X 10-5 

1.0 X 10-^ 5.4 X 10-5 
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Fig. 1. — Spectral maps of the molecular line emission from G324. 201+0. 119. The grid 
spacing is 30" for all lines, except CS(2— >1) which is 45". Offsets are from the reference 
position at q;2ooo = 15^32°^53!2, ^2000 = — 55°56'13". The velocity scale ranges from —120 
to —60 km s~^ . The antenna temperature scale ranges from —0.5 to 3.5 K for CS(2— >1), 
-0.25 to 3.6 K for CS(3^2), -0.2 to 1.0 K for 03^8(3^2), -0.03 to 0.12 K for SiO(2-*l), 
-0.06 to 0.18 K for SiO(3^2), -0.1 to 0.45 K for CH30H(2fe - U), and -0.1 to 0.60 K for 
CH30H(3fe-2fc). 
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Fig. 2. — Spectral maps of the molecular line emission from G328. 307+0. 423. The grid 
spacing is 30" for all lines, except CS(2— >1) which is 45". Offsets are from the reference 
position at ^2000 = 15'^54™06!3, ^2000 = — 53°11'38". The velocity scale ranges from —120 to 
—65 km . The antenna temperature scale ranges from —0.5 to 4.0 K for CS(2^1), —0.5 
to 5.0 K for CS(3^2), -0.2 to 0.9 K for C=^^S(3^2), -0.04 to 0.13 K for SiO(2^1), and 
-0.1 to 0.55 K for CH30H(3fc - 2^). 
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Fig. 3. — Spectral maps of the molecular line emission from G329. 337-1-0. 147. The grid 
spacing is 30" for all lines, except CS(2^1) which is 45". Offsets are from the reference 
position at ctsooo = 16''00™33!4, ^2000 = -52°44'45". The velocity scale ranges from -135 to 
—70 km s^^ . The antenna temperature scale ranges from —0.3 to 4.5 K for CS(2— s>l), —0.6 
to 6.1 K for CS(3^2), -0.2 to 1.4 K for C^^S{3^2), -0.04 to 0.21 K for SiO(2^1), -0.12 
to 0.90 K for CH30H(3fc - 2^), and -0.12 to 0.65 for CH30H(2fe - 1^). 
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Fig. 4. — Spectral maps of the molecular line emission from G330. 949-0. 174. The grid 
spacing is 30" for all lines, except CS(2— >1) which is 45". Offsets are from the reference 
position at 02000 = 16''09™49!1, ^2000 = -51°54'45". The velocity scale ranges from -115 to 
—65 km s~^ . The antenna temperature scale ranges from —0.5 to 4.8 K for CS(2^1), —0.6 
to 6.1 K for CS(3^2), -0.2 to 2.0 K for C''^4S(3^2), -0.07 to 0.75 K for SiO(2^1), -0.15 
to 1.5 K for CH30H(3fe - 2^), and -0.15 to 1.1 for CH30H(2fe - 1^.). 
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Fig. 5. — Spectra observed at the peak position of G324.201+0.119. Transitions are given 
in the upper right corner. The vertical dotted line indicates the systemic velocity of the 
ambient gas. 
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Fig. 6. — Spectra observed at the peak position of G328. 307+0. 423. Transitions are given 
in the upper right corner. The vertical dotted hne indicates the systemic velocity of the 
ambient gas. 
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Fig. 7. — Spectra observed at the peak position of G329. 337+0. 147. Transitions are given 
in the upper right corner. The vertical dotted hne indicates the systemic velocity of the 
ambient gas. 




Fig. 8. — Spectra observed at the peak position of G330. 949-0. 174. Transitions are given 
in the upper right corner. The vertical dotted hne indicates the systemic velocity of the 
ambient gas. 
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Fig. 9. — Contour maps of the CS(2^1) velocity integrated ambient gas emission. Contour 
levels are drawn in steps of 10% and up to 90% of the peak value. The FWHM beam is 
shown in the bottom left corner of the lower left panel. Upper left: G324. 201+0. 119. Peak: 
20.8 K km s'^ . Upper right: G328.307+0.423. Peak: 25.6 K km s"^ . Starting contour 
level: 20%. Lower left: G329.337+0.147. Peak: 28.0 K km s~^ . Starting contour level: 
20%. Lower right: G330.949-0.174. Peak: 39.7 K km s'^ . 
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Fig. 10. — SEST/SIMBA maps of the 1.2-mm dust continuum emission. Contour levels are 
drawn at 5, 10, 20, 30, 40, 50 60, 70, 80 and 90% of the peak flux density. The FWHM 
beam is shown at the bottom left. In all panels, except the upper right, the stars mark the 
positions of the compact radio continuum sources detected by Urquhart et al. (2007). Upper 
left: G324.201+0.119. Peak: 5.28 Jy beam'^ Upper right: G328.307+0.423. Peak: 6.35 Jy 
beam~^. The star marks the position of the bright compact radio continuum source detected 
by Garay et al. (2006). Lower left: G329.337+0.147. Peak: 4.35 Jy beam^^ Lower right: 
G330.949-0.174. Peak: 21.8 Jy beam'^ 
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Fig. 11. — Intensity cuts of the 1.2-mm emission across two orthogonal directions passing 
through the peak position of G324. 201+0.119 (left panels) and G330.949-0.174 (right panels). 
The position angle of the slices are given in the upper left corner of each panel. Dotted lines 
correspond to fits of the observed radial intensity with single power-law intensity profiles. 
The fitted power-law indices are given in the upper right corner of each panel. 
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Fig. 12. — Contour maps of velocity integrated ambient molecular gas emission towards 
G329.337+0.147. Contour levels are drawn at 20, 30, 40, 50 60, 70, 80 and 90% of the peak 
value. The star marks the position of the radio source detected by Walsh et al. (1997). 
Upper: CS(2^1). Peak value: 28.0 K km s"^ . Center: SiO(2-^l). Peak value: 1.45 
K km s-i . Lower: CH3OH 3o 2o A+. Peak value: 3.65 K km s'^ . 



-54- 



10* 



<u 100 



10 



1 1 ' ' 

G324.201 


""1 


1 


= 37 K 
/S= 1.4 




: 


= 5.5 THz 






5^=30" 


/ 

/ 






/ 

/ 

/ 

/ 

/ 

/ 


\ 

\ 

\ 

\ 

■1 


/ 

/ 

: / 
/ 

/ 

/ 




I 

\ 
i 


/ 

1 

■ / 
'■ 1 
1 

1 






1 

1 

1 

1 

' ■ ■ ■ 





1 



10 100 1000 10* 

Frequency (GHz) 



10* 



1000 



-a 

E 100 



1 

G328.307 


l| 1 ' 1 1 1 Mil 

' \ 


= 42 K 


' \ 
' \ 


|0= 1.4 


' \ 
' \ 


i/„ - 5.3 THz 


' \ 
' \ 




' \ 
' \ 


- 


' * 


/ 

/ 

/ 


\ 

\ 


/ 

/ 

/ 




/ 

/ 

/ 

; 

/ 


1 - 

1 

1 


/ 

/ 

/ 

P 

1 




1 
1 

1 


.1 1 . . ' : 



10= 



10 100 1000 10* 105 

Frequency (GHz) 



10* 



G329.337 

Tj = 40 K 
/S= 1.4 
11^ = 5.7 THz 
fl =30" 



CD 

X 
3 

e; 100 



10 



100 1000 10* 

Frequency (GHz) 



1000 



•-3 



d 

<u 100 



X 
3 



10 



1 

G330.949 


1 1 

/ 


: = 37 K 
/S= 1.7 
L/„ - 2.2 THz 

/ 

r / 
/ 

: / 

/ 

/ 

/ 

/ 


1 \ ~ 

/ \ 

1 \ 
/ \ 
1 \ 

\ 

\ 
\ 

\ - 
\ 

\ : 
\ 

\ 


p 

: / 


\ 

\ 

\ 

\ 


/ 

Q 
1 


m\ 


1 

1 

1 

1 

■ 1 

1 

1 

1 


\ 

1 

1 

1 


1 

~ '.1 


1 

1 \'....r, 



10= 



10 100 1000 10* 10= 

Frequency (GHz) 



Fig. 13. — Spectral energy distributions. Stars mark IRAS fluxes, squares MSX fluxes, and 
the circle the SIMBA flux. The long-dashed curve is a flt to the spectrum using two modified 
blackbody functions of the form B^iTd) [l — exp(— (z//z/„)'^)] , with different temperatures. 
The short-dashed line indicates the fit for the colder temperature component (fit parameters 
indicated on the upper left). 
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Fig. 14. — Rotational diagrams for the methanol emission observed toward the massive cores. 
The lines correspond to least squares hnear fits to the observed data. The derived values of 
the rotational temperature and total column density are given in Table [TTJ 
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Fig. 15. — Rotational diagram for the narrow (squares) and broad (stars) SiO emission 
observed towards the G330. 949-0. 174 massive core. The lines correspond to least squares 
linear fits to the observed data. The derived values of the rotational temperature are given 
in the lower left corner. 



